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ABSTRACT 



In this work we give an estimate of the neutrino flux that can be expected from the microquasar Cyg X-3. We calculate 
the muon neutrino flux expected here on Earth as well as the corresponding number of neutrino events in the IceCube 
telescope based on the so-called "hypersoft" X-ray state of Cyg X-3. If the average emission from Cyg X-3 over a 
period of 5 years were as high as during the used X-ray state, a total of 0.1 events should be observed by the full 
IceCube telescope. Using the correlation of AGILE data on the flaring episodes in 2009 Jun.-Jul. to soft X-ray states 
we calculate that the upper limits on the neutrino flux given by IceCube are starting to constrain the hadronic models, 
which have been introduced to interpret the high energy emission detected by AGILE. 

Key words. Neutrinos, Methods: numerical, Cygnus X-3, Microquasars 



1. Introduction 

Microquasars (MQ) are Galactic X-ray bin ary (XRB) 
systems, which e x hibit re l ativist ic radio jets (Mirabel & 



Rodriguez (1999); Fender (2001)). These systems are be 
lieved to consist of a compact object, a neutron star or a 
black hole, and a giant star companion. Mass transfer from 
the giant star to the compact object through the formation 
of an accretion disc and the presence of the jets make them 
similar to small quasars, hence their name "microquasars." 
The analogy may not be only morphological; although there 
is no obvious scaling, it is common thinking that the physi- 
cal processes that govern the formation of the accretion disc 
and the ejection of plasma into the jets are the same for 
both systems. Local, Galactic microquasars may therefore 
be considered as nearby "laboratories," where models of the 
distant, powerful quasars can be tested. The observed ra- 
diation from microquasar jets, typically is in the radio and 
in some cases also in the IR band and it is consistent with 
non-thermal synchrotron radiation emitted by a population 
of relativistic, shock-accelerated electrons. 

The composition of microquasar jets is still an open is- 
sue. The synchrotron emission both in the radio and in 
the IR is consistent with near equipartition between elec- 
trons and magnetic field, which is also implied by min i- 
mum energy considerations (Levinson & Blandford (1996])). 
However, the dominant energy carrier in the jet is presently 
unknown. A possible diagnostic of hadronic jets, namely the 
emission of TeV n eutrinos, has been proposed by |Levinson| 
& Waxman (2001). They have shown that, for typical mi- 



croquasar jet parameters, protons may be accelerated in the 
jet to ~ 10^^ eV, and that the interaction of these protons 
with synchrotron photons emitted by the shock accelerated 
electrons is expected to lead to 1-100 TeV neutrino emis- 
sion. The predicted fluxes are detectable by large, km^-scale 



effective area, high-energy neutrino telescop es, such as the 
o pera ti ng south pole detector IceCube, see Ahrens et al. 
(2003). Distefano et al. (2002) have extended the work of 
jLcvinson &: Waxman, (2001) and predicted the neutrino flux 
for several microquasars whose radio data were available in 
the literature. They found that the largest number of neu- 
trinos from bursting MQs is expected to come from Cyg-X3 
and XTE J1118 +480. 

The AGILE team has made a detailed study of 7- 
ray emission from Cygnus X-3 with the AGILE satellite. 
The AGILE discovery of transient 7-ray emission from 
Cygnus X-3 in 2008 Apr. associated with a specific spec- 
tral state preceding a major radio jet ejection opened 
a new window of investigation of microquasars. Several 



been detected by AGILE and Fermi since 2008 ( 


Tavani 


et al. 


(2009 


); Fermi LAT Collaboration et al. 


(2009 


)). The 



from Cygnus X-3 in the energy range 100 MeV - 50 GeV 
during the periods 2009 Jun-Jul and 2009 Dec.-2010 mid- 
Jun. They found that the 7-ray emission from Cygnus X-3 
is detectable by AG ILE not only during relatively short 



(1-2 day) flares as in Tavani et al. (20091 but also during 
extended perio ds lasting several days or weeks (as during 



2009 Jun-Jul) ( [Bulgarem A. et al.| ( [20T2| )). Detecting con- 
tinuous 7-ray emission during "active phases is of great 
theoretical relevance for the modeling of Cygnus X-3. 

The main mechanism at the origin of the high energy 
emission detected by AGILE can be the Inverse Compton 
scattering of high energy electrons accelerated in the source 
with low energy photon s emitted by synchrotron pro- 



cess (Dubus et al. (2010)). However another possibility, if 
hadrons are present in the jet, is that the 100 MeV - 50 GeV 
emission is due to 7r° decay into TeV 7-rays that trigger a 
cascade to GeV gamma-rays. These 7r° may be produced in 
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the jet by photohadronic interactions or by proton-proton 
collisions between protons in the jet and protons in the 
gaseous surroundings provided by the WR compan ion mass 
outflow ( [Piano et ah] ( |2012| ; [Romero et aT] ( [2003| ) . 

In this paper we use a simplihed jet model to estimate 
the number of neutrinos expected from Cyg X-3 based on 
the X-ray data taken in 2009 during an episode of flares. In 
Section [2[ we describe the simplified model we use for the 
jet of the microquasar and the photohadronic interactions 
happening inside of it. In Section [3[ we derive the expected 
photon and proton densities inside the jet of Cyg X-3, which 
are needed for the simulation of the photohadronic interac- 
tions, from the observed X-ray data. In Section [4[ we then 
use the derived densities to estimate the expected neutrino 
number both in the case that the emission is due to IC and 
in the case that is due to 7r° decay into 7-ray photons. In 
Section |5] we discuss our results. 

2. A simplified jet model 

In this section we want to describe our simplified jet model. 
The model we want to apply to the microquasar Cyg X-3 
was originally used for GRB which are also assumed to orig- 
inate from emitting jets. The main simplification is that we 
totally neglect the dynamic of the propagating jet and as- 
sume that the outflow is static. The high energy protons 
as well as the X-ray and 7-ray photons are assumed to be 
present in one interaction zone and emission will only origi- 
nate from this one emission region. Both the proton as well 
as the photon density are assumed to be isotropic in the 
plasma rest frame of the jet. These assumptions are essen- 
tia l for us to use the phot ohadronic interaction code based 
on Hiimmer et al. (2010b) (model Sim-B), which is an ana- 
lyticai parameterization based on SOPHIA, see |Mucke et al.| 



(2000). Starting from the particle densities the code calcu- 



lates the result of the photohadronic interactions leading 
to secondary particles which subsequently decay into neu- 
trinos. It is based on an analytical parameterization of the 
full photohadronic interaction cross section and includes 
features such as individual treatment of secondary parti- 
cles, helicity dependent mu on decay, as well as n eutrino 

With 



(2010a 



flavor mixing, for details see Hiimmer et al. 
this treatme nt we also incorporate the magnetic field effects 
discussed 



by 

role due to t 



2009), which do play a 



Reynoso fc Romero 

le assumption of neutrino production in the 
internal shocks. 

We assume in concordance with ILevinson fc WaxmanI 



(2001 ) that the photon energy density inside the jet U' can 



be calculated from the emitted luminosity L of the source 
by considering the emission passing the surface of the jet 

L 



47r {P- sm^i9j)) cV^ 



(1) 



with I being the distance of the emission region from the 
central object, V the Doppler factor of the jet, and 9j the 
opening angle of the jet. Note that the boost with the 
Doppler factor is needed since C/^ is the energy density in- 
side the shock frame. Additionally, we know that the photon 
energy density can be defined as 



ds; 



(2) 



with n'^ the photon spectrum in particles per volume and 
per energy. For our photohadronic interaction calculation 



we will need n^. The shape of this spectrum can be ob- 
tained from the photon data while the normalization can 
be obtained using Eq. (d). 

Hence, by knowing the observed (energy integrated) 
photon flux F-y we can calculate the (isotropic equivalent) 
luminosity L = Air Dj^ F-y, and subsequently the photon 
particle density inside the jet. Moreover it is possible to cal- 
culate the energy carried by the magnetic field B as well as 
the energy in protons from ene rgy equipartition arguments . 
Following the considerations of Levinson & Waxman ( 2001 ) 
it is possible to estimate the niagnetic field B' (in G) by 



^Tl — U' 



(3) 



with eb being the fraction of the jets total energy carried by 
magnetic field, Eg the fraction carried by electrons/photons, 
and U'^ from Eq. y\ in ergs~^. Accordingly we can estimate 
that the energy cfensity carried by protons is 



(4) 



with Ep being the fraction of the jets total energy carried 



by protons. The fractions arc estimated to be of the or- 
der of 0.1, see |Distefano et al., ( )2002[) ; 
([2001[). Moreover, we can comiect thi 



Levinson fc Waxman 
e proton spectrum 
to the energy density in protons U'^ by adapting Eq. ^ for 
protons. However, opposite to the photon spectrum, we do 
not have any direct observational information on the shape 
of the proton spectrum. From considerations on the Fermi 
acceleration we know that the spectrum should have a form 
cx E~°''p with ckp ~ 2. Furthermore, we know that particles 
can only be accelerated to finite energies and a cut-off at 
a critical energy Ec should be expectedj^We can estimate 
the critical proton energy by basic considerations on the 
jet properties. In general two limiting scenarios are consid- 
ered: the synchrotron limited case and the escape limited 
case. In the first case the size of the acceleration region is 
sufficiently large to not affect the particle acceleration and 
only synchrotron losses limit the maximal particle energy. 
By comparing the acceleration time to the synchrotron loss 
time it is possible to obtain the critical energy (in plasma 
rest frame) 



E'^^ = 2.01 • 10 



11 



B^ 

Tg 



GeV 



(5) 



with rj the acceleration efficiency and B' the magnetic field 
in the plasma rest frame. For all the calculations in this pa- 
per we have assumed rj = 0.1. In the second case the energy 
of the protons is limited by the maximal energy the par- 
ticles can reach before they escape the source. Effectively 
this is equivalent to comparing the size of the region to the 
Larmor radius of the particles. Hence, this critical energy 
can be estimated by comparing the acceleration time to the 
escape time from the jet which leads to 



E' 



c2 



10% 



10 



14, 



B' 
G 



GeV 



(6) 



The smaller of the two critical energies is the relevant crit- 
ical energy for a given set of parameters. 



The cut-off we assume for our calculations is of the form 
exp (— (_E/i5c)^) , however the actual shape of the cut-off is not 
relevant. Only the position, given by critical energy Ec, is. 
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3. Observational data on Cyg X-3 

The information on Cyg X-3 is comparably well docu- 
mented and there already have been earlier studies on the 
neutri no emission of this source, such as [Distefano et al.| 
(2002). For our study we adopt the following parameters 
for the source: 



Dl 
V 



7.2 kpc 
2.74 , 
12° , 
10* cm 



Di is the distance b etween Ear t h and the source, with 
the result taken from Ling et al. (20091. The value of the 
Doppler factor T) is the result of an assumed Lorentz factor 
r = 1 . 70 and a viewing angle = 14°, tak en from |Distefano| 
et al.| (|2002[ ) ; |Mioduszews ki et al .| (|200lD. The jet op ening 
angle^j is also taken from Mioduszewski et al. ( 2001| ). The 
value for the radius of the emission region I is set to the as- 
sume d collision radius of internal shocks in a microquasar, 
as in Levinson & Waxman (2001). 

In a second step we need to take into account the kine- 
matics of the photohadronic interactions to estimate which 
energy range of photons can act as the target photons 
for our highest energy protons. Assuming a maximal pro- 
ton energy of 10* GeV we can estimate that X-ray pho- 
tons with several keV energy are needed for producing the 
A-resonance. F o r this reason we will use X-ray data from 
Koljonen et al. (2010) as our target photons. Since we are 
interested in the phase directly before a flare we will use 
data for the hypersoft state from said paper. The used data 
ranges from 3.5 up to 102 keV, with the integrated luminos- 
ity we obtain being = 3.77- 10'^'' erg s^^. With help of 
Eqs. ([T]) and ^ we can normalize our target photon density. 
Furthermore we can estimate that for the given the re- 
sulting magnetic field from Eq. ([s]) is -B' = 8.8 • 10'"' G, while 
the resulting maximal proton energy is E'^^^ — 5.5 • 10'' GeV, 
given by Eq. (6|. Hence we can normalize the proton den- 
sity as well through Eq. ([4]). The density, which we ob- 
tain using this method, corresponds to an injection rate of 
3.3 • 10^^ protons cm~^ s~^. 



4. Expected neutrino flux 

With the normalized spectra we obtained in the previous 
section we can calculate the shape of the neutrino flux 
at Earth after flavor mixing. Especially the shape of the 
muon neutrino flux is relevant for the current neutrino ex- 
periments such as IceCube. Our photohadronic interaction 
cross section includes several contributions apart from the 
A-resonance such as higher resonances, direct production 
of pious (i-channel), and high energy processes leading to 
multiple pious. For the flavor mixing we still assume a sce- 
nario with 9i3 = eve n though Daya B ay, see |An et alT] 
(2012|, and RENO, see |Ahn et all ( |2012h , have ruled this 
out. However, it was already shown that the effect of differ- 
ent values of the neutrino mixing angles inside the current 
uncertaint y on the neutrino flux prediction is comparably 
small, seelBaerwald et al. (2011), and can be neglected at 



for GRB since we apply the same assumption of a neu- 
trino flux originating from internal shocks. Hence we can 
see the same splitting into a double peak from muon and 
pion decay plus an additiona l component from kaon decay. 



as m [Baerwald et al. (2011). It is still significantly below 
the detection limit of IC59 during the time the flares were 
recorded (assumed exposure of 61 days, black solid curve), 
and if one extrapolates to the full IceCube detector and 5 
years of exposure the neutrino flux would just start to be 
in the detectable range (black dashed curve). Note that the 
here depicted limits are based on the solid-angle-averaged 
effective areas at final cu t level of the time- inte grated point- 
like source search from Abbasi et al. (2011). We use the 
band for S = (30°, 60°) from the left plot of Fig. 8 of 
said reference, which should be is suitable for Cyg X-3 
{S = 40° 57'). Even though this effective area is already 
at the final cut level it still needs to be considered that the 
neutrino signal needs to be distinguished from the back- 
ground of atmospheric neutrinos. There are however possi- 
bilities to improve the cuts, such as ad ditional cuts in the 
timing as done in Abbasi et al. (2012) for Cyg X-3 flares. 
Additionally we have scaled the effective area up from the 
given values for the 40-string conflguration to the 59-string 
(multiplied by 1.3) and 86-string configuration (multiplied 
by 2.4). Therefore the actual effective areas for IC59 and 
IC86 point source searches may actually be slightly differ- 
ent, but were not publicly available during our work. 

Moreover, it can be seen that the predicted neutrino 
flux above 10^ GeV is already one order lower than the 
peak value and the flux drops off even more above PeV 
energies. Therefore, it should be unlikely that Cyg X-3 is a 
possible source for the t wo uniden t ified P eV cascade events 
in IceCube, see talk by Ishihara (2012) at Neutrino2012, 
as events at lower energies, i.e. between 10^ and 10^ GeV, 
are far more likely. Furthermore, the number of events can 
be calculated from the flux prediction and the detector pa- 
rameters using the formula 



N = J dEAf{E) texp 



d N{E) 
dE 



(7) 



with Alf{E) being the energy dependent effective area (of 
IceCube), tcxp the exposure, and dN{E)/dE the (muon) 
neutrino spectrum (on Earth, after flavor mixing). We ob- 
tain that we would expect 0.002 events in IceCube for the 
time of the flares in June and July of 2009 (59 strings) 
and a total of 0.12 events over 5 years in the full detec- 
tor. Especially this flrst result of no events in IC 59 is 
(in a sense) reassuring as it is consistent with the current 
IceCube data, which so far does not suggest any UHE neu- 
trinos apart from the atmospheric ones, see Abbasi et al 
(2012). Moreover, it should be possible to constrain ele- 



ments of this simple microquasar model with in 5 to 10 
years of data taking with the full detector if the flares con- 
tinue. We are aware that this model is very simple and does 
not even try to explain the dynamic of a microquasar jet. 
However, as neutrino measurements are at this stage not 
even close when it comes to the time resolution compared 
to photon measurements, it is reasonable to average out 
most of the variations. 



this stage. In Fig we plOTTe expected muon neutrino ^^^t^;^^, ^ 



(and antineutrino) flux (in GeVcm^ s^ ) for Cyg X-3 in 
this simple model of a microquasar jet as a blue solid curve. 
The shape of the spectrum is similar to the ones predicted 



We now want to test how many neutrinos would be ex- 
pected if the observed 7-ray emission by AGILE was actu- 
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6 10" 




^JC86 (5 years) 




Fig. 1. In this plot we show the expected neutrino flux from 
Cyg X-3 derived for a simplified model of the microquasar 
j et, ba sed on hypersoft X-ray data from |Koljonen et al. 
(2010) (blue solid curve). For comparison we give the (ex- 
trapolated) differential flux limits for IC59 during the 61 
days of the flaring episodes (solid black curve) as well as an 
extrapolation of the expected one for IC86 (full detector) 
after 5 years of data taking. As one can see both do not 
reach the flux of this single microquasar, even though the 
full detector could with a longer exposure. The differential 
limits are based on the solid-angle-averaged effective areas 
at final cut level of the time-in tegrated point-like source 
search from Abbasi et al. (2011 ). 



ally coming from the decay of photohadronically produced 
TT*^ into photons. The photons from such decays would have 
to cascade down to lower energies and may lose a part of 
their energy during this process. However, here we only 
want to estimate the number of expected neutrinos for the 
case that the observed GeV-emission is due to 7r° decays. 
Hence, it is sufficient to assume energy conservation for 
the calculation. With this approach it is possible to set a 
lower bound on the expected amount of neutrinos. In case 
of energy losses during/after the cascading process an even 
higher amount of original tt'^ (and consequently 7r+) would 
be needed, and hence an even higher neutrino flux would 
be expected. If one assumes that 



IOGgV 



dE-y E-y 



AN^{E^, 



10"' GoV 



dK 



d-E^ E^ 



AE^ 



(8) 



lOMoV 



OGcV 



with '^^^^^'^ being the observed photon spectrum from 



Bulgarelli, A. et al. 



(2012) while 



is the rcsult- 

mg spectrum ot high energy photons before cascading to 
lower energies. The high energy spectrum has been calcu- 
lated numerically with the photohadronic interaction code, 
which also has been used for calculating the neutrino spec- 
tra. With this method we obtain that the amount of pro- 
duced tt" needed to match the observed 7-ray emission is 
about 2400 times larger than our prediction for tt'^ from the 
calculation in the previous section. As a consequence the 
nominal amount of expected neutrino events would reach 



about 5.2 events for the 61 days of flaring in 2009. This 
would still need to be distinguished from the O(IOO) back- 
ground eventf]^ If the used cuts achieve such a high preci- 
sion, then it should be possible to rule out if the observed 
flares are due to the decay of tt'^. 

On a similar note it can be stated that models based 
on pp-interactions might already be more strongly con- 
strained than the photohadronic interaction models. Even 
though we cannot give a detailed estimate on the level of the 
proton-photon collisions for the proton-proton interactions 
we can still give a zeroth order estimate. It is commonly 
assumed that pp-interactions give a neutrino result which 
is 4/3 higher than for p^. Hence this very rough estimate 
gives that 0.16 neutrino events should be expected over 
5 years or 6.9 events during the flaring episodes for 7-ray 
flares from Tr^-decay. A more detailed calculation of the neu- 
trino flux prediction from proton-proton interactions, with 
all features included for the photohadronic interactions, 
is currently beyond the scope of this paper e ven though 
work in this dir ection was already done, see e.g. |Reynoso fc| 
Romero (2009). A prediction based on a detailed particle 
physics treatment of pp-intcractions should be a goal for 
the future. 



5. Discussion 

In this paper we have investigated the possibility that nii- 
croquasars may be the sources of high energy neutrinos. 
In particular we have estimated the neutrino flux expected 
from Cygnus X-3. Starting from the X- ray data for the hy- 



persoft state from 
back to the partic 



Koljonen et al. (2010) we have calculated 
e densities of protons and photons inside 



the jet of Cyg X-3 in the assumption of a simple geometrical 
model for the jet. The reason for choosing this data set was 
to use a data for a phase bef ore a radio flare, as suggested 
by lBulgarelh, A. et al ||2012[ ). This state also fulfills the ob- 
served correlation of AGILE 'j-ray flares in June- July 2009, 
and November 2009- July 2010 with soft X-ray states and 
episodes of decre asing or non-detectable har d X-ray emis- 
sion reported by Bulgarelli, A. et al. (2012). However, in 
principle any other state could also be used tor the calcula- 
tions, as long as we are able to derive the particle densities 
inside the jet. We then used the densities to compute the 
expected neutrino emission from the jet using a numerical 
code which incorporates the full photohadronic interaction 
cross section, individual treatment of secondary particles 
(including losses), and flavor mixing of the neutrinos. The 
expected muon neutrino flux was then compared to the sen- 
sitivity of IceCube (59 strings) during the 61 days of flaring. 
The expected number of 0.002 events is in concordance with 
the non-detection of any neutrinos from Cyg X-3 during 
that period. Assuming an extended period of flaring and a 
full IC86 detector we expect about 0.1 associated neutrino 
events in 5 years of data taking. Moreover, the shape of 
the neutrino flux disfavors a detection at PeV before seeing 
events at about lOTeV. Nevertheless, these predictions are 
still subject to quite large uncertainties due to the error on 
the jet parameters as well as the used extrapolated effec- 
tive area of IceCube. Therefore even slightly lower amounts 
of observed neutrinos may not directly contradict the ba- 



^ The amount of events is approximated from scal ing the 
amount of e vents at final selection level from Fig. 4 of | Abbasi] 



et al. (20111 down to 61 days (from 375.5). 
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sic model of microquasars accelerating protons. Still, in the 
best case it should be possible to see some events from 
Cyg X-3 with several years of data taking with the full 
IceCube detector. Moreover, we also tested the hypothesis 
that the observed 7-ray emission is due to the decay of 7r° 
from photohadronic interactions into photons. For this we 
compared the integrated energies in the observed photons 
detected by AGILE to the energy in photons from 7r° de- 
cay. In the A-resonance approximation the production of 
77° is directly connected to the production of tt"*", and this 
concept in principle does not change even for a more de- 
tailed particle physics treatment of the interactions. Using 
these calculations we obtained that the amount of energy in 
7r° needed to explain the AGILE observations would have 
to have been so high that the number of expected neu- 
trino events would reach about 5 events during the 61 days 
of flaring. There have, however, been no reported neutrino 
events in Ice Cube which could be associated with Cyg X- 
3 so far, see lAbbasi et al. (2012). We can therefore start 



to rule out that the observed 7-ray emission is due to the 
decay of ir^ from photohadronic interactions by combin- 
ing the photon and neutrino data in the coming months. 
Especially the point source analysis of IC59 would be of 
great interest in this regard. Moreover, based on the pho- 
tohadronic estimates and the estimate that the result for 
proton-proton collisions is larger by a factor of 4/3, it can be 
assumed that pp-interactions are even stronger constrained 
than j>7-interactions. However, a real simulation of the neu- 
trino flux prediction from pp, based on a detailed particle 
physics treatment, is still needed for a conclusive picture. 
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